Optical response of many-polaron systems 
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Abstract 

Exact results for the density of states and the ac conductivity of the spinless Holstein model at finite carrier density are 
obtained combining Lanczos and kernel polynomial methods. 
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Optical measurements have proved the importance of 
electron-phonon (EP) coupling and even polaron effects 
in several important classes of materials, including one- 
dimensional (Id) MX chains, quasi-2d cuprate superconduc- 
tors, and 3d colossal magnetoresistance manganites. In all 
these materials, a noticeable density of (polaronic) charge 
carriers is observed, which puts the applicability of single- 
polaron theories into question, particularly in the region of 
intermediate EP coupling strength and phonon frequency. 

Recently, the photoemission spectra of many-polaron sys- 
tems have been investigated in the framework of the Id 
spinless Holstein model, 
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describing tight-binding (f) electrons coupled locally (g) to 
Einstein phonons (wo), where c] (b\) denote the correspond- 
ing fermionic (bosonic) creation operators, and n, = c? i c i . 
Most notably, provided that the EP coupling is not too strong, 
a density driven crossover from large polarons to weakly 
dressed electrons has been found to occur [1], In the mean- 
time, this result was corroborated by cluster perturbation 
theory (CPT) [2]. 
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Here, we use Lanczos diagonalization and kernel poly- 
nomial expansion methods [3] to study the (linear) optical 
response of a many-polaron system to an external (longi- 
tudinal) electric field, Rea(uj) = D8(u) + a rcs (uj), the 
regular part of which is given by 
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where dJ mn = E m — E n , E n (|n)) are the eigenvalues 
(eigenstates) of our iV-site coupled EP system with at most 
M phonons, and j — i et 2i( c I c i+i — c l+i c i )■ in addition, 
we calculate the partial densities of states (DOS) p + {u) and 
p~(u), which are obtained from the fc-integrated single- 
particle spectral functions 
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c t = C L c k ~ c fe> wnere (u>) [A^(ui)] is related to the 
[inverse (I)] photoemission (PE) of an electron. 

Figure 1 displays selected numerical results for the optical 
properties of the Id model (1) for g 2 = 5, wo/i = 0.4 and 
various characteristic particle densities. 

Starting with a single electron (n — 0.1), we notice from 
the DOS that there is a polaron feature at the Fermi level Ef 
(cf. the jump in the integrated DOS at Ef and the negligible 
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Fig. 1. (Color online) Left panel: Single-particle DOS p~ (cu) (red dashed lines) and p 4 " (a?) (black solid lines) for g 2 — 5, U)o/t — 0.4, and particle 
densities n — 0.1 (one electron), n — 0.3, and n — 0.5 (half filling). Results are for a ten-site system with periodic boundary conditions and M < 15 
dynamical phonons. The homogeneous q — lattice displacement was treated separately [5], leading to truncation errors < 10 -4 . Dashed (magenta) 
lines give the integrated DOS. CPT data for p^fu) included for comparison (blue lines) were partly extracted from [2]. Right panel: Regular part of the 
optical conductivity (7 Te£ (uj). (Green) dashed-dotted lines: analytical strong-coupling result a TCR (cu) — oq n (uo^) -1 uj~ x cxp[(a; — 2g 2 loq) / (2guJo)] 2 
(<to = 8) [4]. 



spectral weight of p~ (cu < Ef)). However, in view of the in- 
termediate EP coupling strength chosen, A = g 2 u>o/2t — 1, 
the polaron is rather extended (large). Consequently, the po- 
laron band is not far separated from incoherent excitations, 
and cr rGg (a>) strongly deviates from the analytical strong- 
coupling result. In particular, the maximum in <j rcg occurs 
well below the small polaron value 4A. 

At finite carrier density, the system shows diffusive trans- 
port. The polarons are dissociated and the remaining elec- 
tronic quasiparticles are scattered by virtual phonons. The 
peaks in the PE part of the spectrum, which for n = 0.1 
had reflected the Poisson-like distribution of phonons in 
the ground state, now broaden significantly, and ultimately 
merge with the IPE part into a wide incoherent band (see 
also the continuous increase of the integrated DOS) [1]. 

In the half-filled band case, the model has a symmetry- 
broken insulating charge density wave ground state accom- 
panied by a distortion of the lattice. This is because the 
EP coupling exceeds the critical interaction strength for the 
Peierls instability: g > g c (ujo/t = 0.4) ~ 1.9 [5]. Accord- 
ingly, we find a gap feature in the DOS and a clear optical 
absorption threshold. Note that the CPT does not reproduce 
the very small gap in the DOS, as it "interpolates" between 
the wavevectors of the finite cluster. 

In conclusion, increasing the particle density in the 



Id spinless Holstein model at intermediate EP coupling 
strengths, we observe a crossover from a polaronic system 
to a metal composed of weakly dressed electronic quasi- 
particles and finally to a Peierls insulator at half filling. 
These transitions are reflected in significant changes of the 
optical spectra. 
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